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Abstract. In "Decoherence of macroscopic closed systems within Newtonian 
quantum gravity" (Kay B S 1998 Class. Quantum Grav. 15 L89-L98) it was argued 
that, given a many-body Schrodinger wave function . . . ,xjv) for the centre-of- 

mass degrees of freedom of a closed system of N identical uniform-mass balls of mass 
M and radius R, taking account of quantum gravitational effects and then tracing 
over the gravitational field amounts to multiplying the position-space density matrix 
p( x i, . . . , xjyjxi, • • • jXjy) = V'( x i7 ■ ■ ■ i x Af)V'*( x 'i: ■ ■ ■ j x at) by a multiplicative factor, 
which, if the positions {xi, . . . ,xjy;x^, . . . ,n' N } are all much further away from one 
another than R, is well-approximated by 

(ru ( J =s =ki ) il<j ( \^-ZtV-^ \)y 24M ■ Here we show that if each uniW 

mass ball is replaced by a grainy ball or more general-shaped lump of similar size 
consisting of a number, n, of well-spaced small balls of mass m and radius r and, in 
the above formula, R is replaced by r, M by m and the products are taken over all Nn 
positions of all the small balls, then the result is well-approximated by replacing R in 
the original formula by a new value R e Q- This suggests that the original formula will 
apply in general to physically realistic lumps - be they macroscopic lumps of ordinary 
matter with the grains atomic nuclei etc. or be they atomic nuclei themselves with 
their own (quantum) grainy substructure - provided R is chosen suitably. In the case 
of a cubical lump consisting of n = (21 + 1) 3 small balls [l > 1) of radius r with centres 
at the vertices of a cubic lattice of spacing a (assumed to be very much bigger than 
2r) and side 2£a we establish the bound e _1/3 (r/a) 1 /"ia < i? cff < 2 v / 3(r/a) 1/ "£a. 
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1. Introduction and Results 

In P , one of us argued that several of the puzzles inherent in our current understanding 
of quantum and quantum-gravitational physics would appear to find a natural resolution 
if one were to postulate that quantum gravity is a quantum theory of a conventional 
type with a total Hilbert space -f^otal which arises as a tensor product 

-^total = -^matter ® -^gravity 

of a matter and a gravity Hilbert space, and with a total time-evolution which is unitary, 
but that, while, as would usually be assumed in a standard quantum theory, one still 
assumes that there is an ever-pure time-evolving "underlying" state, modelled by a 
density operator of form 

Ptotal =!*><*!. 

^ G -f^totab a ^ eacn " ms tant of time", one adds the new assumption that the physically 
relevant density operator is not this underlying density operator, but rather its partial 
trace, p ma tter' over -^gravity With these assumptions, an initial underlying state of a 
closed quantum gravitational system with a low degree of matter-gravity entanglement 
would be expected to become more and more entangled as time increases - equivalently 
Pmatter will De subject to an ever-increasing amount of decoherence. In fact the von- 
Neumann entropy of p ma tter would be expected to increase monotonically, thus (when 
the theory is applied to a model for the universe as a whole) offering an explanation 
for the Second Law of Thermodynamics and (when the theory is applied to a model 
closed system consisting of a black hole sitting in an otherwise empty universe) offering 
a resolution to the Information Loss Puzzle. (It offers both of these things in that, 
by defining the physical entropy of P^otal ^° ^ e ^ ne V °n-Neumann entropy of p ma tter' 
one reconciles an underlying unitary time-evolution on -f^otal with an entropy which 
varies/increases in time.) 

A second paper, [2] , by the same one of us investigated the implications of the theory 
proposed in pQ for the decoherence of ordinary matter in the non-relativistic, weak- 
gravitational-field regime. In particular, [2 addressed the question of the appropriate 
description of the state, at some instant of time, of the center-of-mass degrees of freedom 
of a system of N bodies, which, for simplicity, were taken to be identical. In ordinary 
quantum mechanics, such a state would be described by a many-body Schrodinger wave 
function ^(x 1 , . . . , x N ) in the usual matter Hilbert space -ff ma tter' consisting of the 
(appropriately symmetrized, according to whether the bodies are treated as fermions or 
bosons or neither) iV-fold tensor product of L 2 (R 3 ). [2j argues that such a ip needs to be 
replaced by a total (entangled) matter-gravity state vector in a total matter-gravity 
Hilbert space which is the tensor product of -^matter (defined as above) with a suitable 
Hilbert space -^gx- a vity ^ or the modes of the quantized linearized gravitational field. 
Specifically, thinking of this tensor-product Hilbert space as the set of (appropriately 
symmetrized) square-integrable functions from R 3N to -^g rav ity 1^0 is taken, in 0, to 
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be the function 

(x 1 ,...,x iV )^(x\...x N )| 7 (x 1 ,...x N )> 

where ^(x 1 , . . .x N )) is a certain (non-radiative) quantum coherent state (introduced in 
PI) of the linearised gravitational field describing the Newtonian gravitational field due 
to the simultaneous presence of a body centred on each of the locations x 1 ,...x N . 
(We remark that, in |2j, for the purposes of calculating these coherent states, the 
bodies are modelled as classical mass-distributions.) In the non-relativistic, weak- 
gravitational-field regime, the projector p^otal = on ^° this \^f), is thus taken 
to be, to a very good approximation, the correct description of the "underlying" state 
of quantum gravity describing our many-body system in the sense described in the 
previous paragraph. The physically relevant density operator is therefore given, in this 

approximation, by the partial trace p ma tter °f this Ptotal over -^gravity wmcri ; in 
position space is clearly given by 

Pmatter( x i> • • • > x n] x d • • • , x a0 

= V(xi, . . . ,x Ar )^*(x' 1 , . . . ,x^)A^(xi, . . .jX^xi, . . . ,x^) 

where the multiplicative factor .M(xi, . . . , xtv; x' 1; . . . , x^) is defined by 

M(xi, ■ • .,xjv;xi, ...,x^) = (7(x 1; . . . ,x 7V )|7(x' 1 , . . .,x.' N )). 

(In [2 , M. was written e~ D where D was called the decoherence exponent.) Moreover, 
it was shown in [2] that, in the case the bodies are taken to be balls with constant mass 
density and (we work in Planck units where G — c — H — 1) total mass M and radius R 
(i.e. in the case the classical mass distributions representing the bodies are taken to be 
such balls) and in the case all the positions {xi, . . . , xjv; x' 1; . . . , x^} of their centres of 
mass are much further away from one another than R (we shall call this the well-spaced 
regime) J, the multiplicative factor is well approximated by M. a where M. a is given by 
the explicit formula 

N N / 1 ' _ II _ / |\ -12M 2 

Ai4x,.,x,x;,..,x; ) = nn ! */ (D 

/=1J=1 VI X / x Jll x / x jI/ 

where it is to be understood that, in the cases I = J, the terms in the denominator 
|x/ — xjIIx'j — Xj| are to be replaced by R 2 . 

To summarize, taking into account gravitational effects and then tracing over the 
gravitational field, as the general theory of £Q prescribes that we should do, has, in the 
non-relativistic weak-gravitational-field regime, and on the assumption that the bodies 
are uniform mass balls, an overall effect which is equivalent to multiplying the usual 
quantum mechanical (pure) position-space density matrix ^(xi, . . . , xat)?/)*(x / 1 , . . . , x^) 

| In many applications, one expects the region of (configuration space) x (configuration space) where 
this condition doesn't hold to be so small in comparison to size of the region where ipip* is 
significantly big that no significant correction would be needed to results calculated on the assumption 
that the multiplicative factor is well-approximated by the M. a in all of (configuration space) x 
(configuration space). 
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by a multiplicative factor, Ai, which, in the well-spaced regime, is well-approximated by 
Ai a given by the formula (0) and it is the product of this pure density matrix with the 
multiplicative factor which is to be regarded as the physically relevant density operator. 
We remark that the formula (0) can also be written 

(. . \ -2AM 2 
n ( |xy ~/ k1 ) n ( 
K ^ ' KJ ^ 

where the first product is taken over all K from 1 to N and the second product is taken 
over all I and J from 1 to N, which satisfy I < J. We also remark that, in the case 
of a wave function for a single ball state, this prescription amounts to multiplying the 
density operator ■0(x)^(x / ) by the multiplicative factor 

^ a (x;x') = (|x-x'|/i?)- 24M2 . 

The problems we wish to address in the present paper are that, in the derivation, [2j, 
of the bodies are assumed to be constant mass-density balls whereas in applications, 
the bodies one is interested in will actually typically be "grainy" , and also not necessarily 
spherical. In fact, in one application (cf. the discussion of the Schrodinger-cat-like states 
in of the above formulae, the balls are interpreted as macroscopic balls of ordinary 
matter and real ordinary matter is of course grainy in that it is made out of atoms 
etc. and we may also be interested in lumps of ordinary matter with shapes other than 
balls. In another application, the above formulae (or rather their obvious generalization 
to states of many bodies where the masses and radii are not all equal) are interpreted 
as telling us how the standard non-relativistic quantum mechanics of closed systems of 
large numbers of atomic nuclei and electrons etc. gets modified according to the theory 
developed in pQ and [2|. In this latter application, the balls are taken to be models of 
nuclei and electrons etc. and again, of course, real nuclei and electrons are not actually 
uniform density balls, but will have their own grainy substructure and also need not be 
spherical. 

However one might hope that if const ant- density balls were replaced by grainy, and 
possibly non-spherical, (classical) "lumps", then the formula for the multiplicative 
factor might nevertheless remain approximately correct provided we replace R in (PJ by 
a suitable effective radius R e Q- One might further hope that, in the case of a grainy 
ball, i? e g would turn out to be of the same order of magnitude as R and in the case 
of a grainy lump of some other shape, to be of the same order of magnitude as some 
measure of the lump's typical linear size. (We shall refer below to the lump's diameter 
without attempting to give a precise general definition to this notion.) If these hopes 
were fulfilled, then one might say that the formula for M. a is robust under (classical) 
changes in graininess and in shape. 

We remark that we are continuing to assume here that our grainy, non-spherical, 
situations can, as far as their gravitational effects are concerned, still be modelled within 
the formalism of [2] , as classical (albeit now no longer uniform and no longer spherical) 
mass distributions and we admit that, in principle, we should of course presumably work 
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within an extension of the formalism of 2 which allows for truly quantum (and also 
relativistic) descriptions of graininess (and non-sphericity). (Especially, such a quantum 
relativistic description may well be important in the second of our applications, to the 
grainy substructure of the proton etc.) While we shall not attempt to explore such 
quantum notions of graininess and non-sphericity in this paper, if robustness holds in 
the above sense under classical changes in graininess and in shape, it might also seem 
reasonable to guess that our formulae will remain robust even if graininess and changes 
in shape were taken into account in such a correct quantum way. 

In any case, the purpose of the present paper is to investigate whether and/or 
under what circumstances, these hopes are fulfilled in the case of one specific class of 
classical grainy models. Namely, models in which every one of the balls of mass M 
and radius R in formula is replaced by a collection of n little balls, each of mass m 
and radius r fixed at definite positions inside a lump-shaped region (i.e. some smooth 
closed bounded region of R 3 ) of diameter approximately equal to 2R so that its centre 
of mass is located where the center (i.e. one of the positions xi, . . . , xjv) of the ball it 
replaces was located. We shall assume that the configurations of the balls within each 
of the lumps are Euclidean-congruent to one-another. However we shall not restrict the 
relative orientations of the different lumps; they are allowed to be arbitrarily rotated 
with respect to one another. 

We remark that obviously we have nm = M. Here we have in mind choices of 
values for the various parameters such that 

r <C any inter-small-ball spacing and R <C any inter-lump spacing, (2) 

where, by "any inter-lump spacing" we mean any of the distances between pairs of 
distinct elements of the set {xi, . . . , x^; x' 1; . . . , x.' N } and (assuming from now on that 
a definite system has been adopted for enumerating the little balls - as Xu, . . . ,xi n 
etc. - within the Jth lump etc.) by "any inter-small-ball spacing" we mean any of the 
distances between pairs of distinct elements of the set {x/i, • • • , Xj n ; x' n , . . . , x^ n } for 
some/any I. Then (in all but an unimportantly small volume of (configuration space) x 
(configuration space) - cf. footnote above) formula (0) will clearly get replaced by 

N n N n / 1 / II _ / I \ ~ 12m2 

. , x n; x ;, , = n n n n bi^*":^ < 3 > 

i=i i= i J=1 j= i V |X/J Xjj 1 |X/i Xj j 1 / 

where now it is to be understood that in the cases where I = J and i = j, the terms in 
the denominator — \x' u — XjJ are to be replaced by r 2 . 
In line with (J2J), if we make the replacements: 

— xjj — > Xj — xj, and x/j — xjj — > x/ — xj except when I = J 

(and similarly with primed and unprimed quantities interchanged) in (JHJ), then we clearly 
expect to get a good approximation§ to A^f ram ^. Making these replacements, we 

§ In fact this will clearly be true in the sense that if we scale the positions of the centres of mass of 

the lumps while not scaling the lumps themselves, then A^f iam ^/A^ a tends to 1 as the scale tends to 
infinity. 
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immediately get 

^f ain y(x 1 ,...,x JV ;x' 1 ,...,x' JV ) 



ILIL, 


l X 7 - X J 


|„ 2 


|x 7 - jdj 


|n 2 




X/ - Xj| 


n 2 | 




n 2 



n 



1,3 



Xj - X 



2N \ -12m J 



mil, 


|Xj - Xjj 


l x / - x 'jl 




X/ - Xj 


X ',- X SI 



n 



where the final product involves all the positions Xj of little balls inside any single lump 
- by our assumption of Euclidean congruence, it doesn't matter which one - and ranges 
over all values of i and j from 1 to n except that it is to be understood that, in this final 
product, in the cases i = j the denominator is to be replaced by r 2 . Clearly another 
way of saying exactly the same thing as this is that TVf S rain y j g approximately equal to 
-M a , as given by the formula (JI} except that the proviso that, in the cases I = J, the 
terms in the denominator |xj — xj| |x'j — x'j| are to be replaced by R 2 should be replaced 
by the proviso that, in the cases I = J, the terms in the denominator |xj — xj| |xj — x'j| 
are to be replaced by where i? e g is defined by 



%= nn 



X^ Xj I 



J 



where the products each go from 1 to n and it is to be understood that, in cases where 
i = j, the terms |xj — Xj| are to be replaced by r. Equivalently, we may write 

(where the product is now over all i and j from 1 to n satisfying i < j). Thus the first 
part of our hope (i.e. that replacing our uniform density balls by our grainy lumps can 
be well-approximated by replacing R by a suitable R g q) is fulfilled. 

We now turn to discussing the second part of our hope, namely whether, and/or 
under what circumstances, R q q will turn out to be of the same "order of magnitude" 
as the diameter of our lump. We have not been able to answer this question in general 
but, instead, content ourselves with analysing one particularly simple special case of our 
model. Namely, where each lump is a cube consisting of n = (2£ + l) 3 small balls of 
radius r with centres at the vertices of a cubic lattice of spacing a and side 2£a. To spell 
out precisely what we mean and at the same time set up a useful notation, we assume 
the ball-centres of any one of these cubes to be coordinatizable so that they lie at the 
positions (iia, i 2 a, i 3 a) where ii, i 2 , and i 3 are integers which each range between 1 and 
2£ + 1. Clearly the total number of balls, n, in the cube will be (2£ + l) 3 and we note 
that there will be a ball at the centre (with coordinates ((£ + l)a, (£+ l)a, (£+ l)a)). In 
line with (J2J), we assume 

r a and R -C any inter-lump spacing. 



The robustness of a many-body decoherence formula of Kay 



7 



(Of course, we require r < a/2 for the balls to fit into the lattice at all!) 
For such a lump, Q can clearly be rewritten in the form 



i 

r \ V* / -n- X,; — Xj 



n 



(5) 

u / 

V t<J / 

(here we continue to assume some system has been adopted for numbering our balls 
from 1 to n — (2£ + l) 3 ) and we notice that, for many values of the pair (r/a,n), the 
prefactor, (r/a) 1 /" will itself be of order 1. This will be the case if, as is relevant to a 
ball or lump of "ordinary matter", we take r/a to be of the order of the ratio of the 
radius of the proton to the Bohr radius (i.e. around 10~ 5 ) even for £ = 1 (n = 27) 
and in fact, if £ is 2 (n = 125) or more (and it will be much more if we are taking our 
cube to be a model for a macroscopic lump of ordinary matter as considered in [2]) then 
this prefactor will in fact be a number close in value to 1. On the assumption that this 
prefactor is of order 1 (or very close to 1) our question then reduces to the question 
whether the quantity 



n 



is of the same order as £. Choosing units such that a = 1, this amounts to a question 
about the quantity || 




q(£)= n^-^i ) ( g ) 

\i<j 

where the product is over all pairs of distinct points of the cubic lattice of points with 
integer coordinates (11,12,13) where i\, %i and 23 range between 1 and 2£ + 1 (from now 
on we shall often simply call such a finite cubic lattice of points a "cube") and, again, 
we assume that some system for numbering these points from 1 to n — (2£ + l) 3 has 
been adopted. We formulate this question in a precise way by asking whether one can 
find numbers, Ci, c 2 , of order 1 such that 

at <Q{£)< c 2 L (7) 

It is easy to see that one can satisfy the upper bound by choosing c 2 = 2 a/3. To see 
this, notice that each term in the product in © is obviously less than or equal to 2\f?>£ 
since this is the length of the body-diagonal of the cubic lattice and moreover there are 
n(n — l)/2 terms in the product. We thus have that Q(£) < (2y/3£) n( - n ~ 1 ^ n2 which, since 
£ > 1, is clearly less than 2^/3£. 

However, as far as we can see, to show that there exists a (non-zero) positive C\ such 
that c\£ is a lower bound is not completely trivial because Q(£) is the 2/n 2 power (i.e. 
2/{2£ J r l) 3 power) of a product of numbers which range in magnitude from 1 to numbers 
of order £ (the largest of the numbers will of course be 2y/3£). Nor, to our knowledge, 

|| It is necessary to bear in mind, here and elsewhere, that n is shorthand for (2£ + 1) . 
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does it easily follow from any existing standard mathematical result^ However, we have 
succeeded in finding the (albeit probably not best-possible) number c\ = e -1 / 3 for which 
we can prove that it does hold: 

Proposition. e~ 1/3 £ < Q{£) < 2\fl£. 

(Actually, one can show H , by similar methods to those in the proof of the lower bound 
below, that if we omit £ = 1 the upper bound of 2y/3£ can be reduced to 2.611.) We 
present our proof of the (lower bound in) the proposition in a separate section below. 
Putting this proposition together with (|5|). we thus have the following theorem. 

Theorem. For a lump consisting of n = {21 + l) 3 small balls of mass m and radius 
r, centred at the vertices of a cubical region (as specified above) of side 2£a of a cubic 
lattice of spacing a, 

e- x l\T/a) x f n la < R eff < 2^{r / a) l,n £a. 

(Here, we recall that a is of course assumed to be greater than 2r and R e jj, when inserted 
in place of R in (0), is expected to give a good approximation to the multiplicative factor 
M grainy o/ p z/a>2r .j 

We remark that it is easy to calculate Q(£) numerically for small values of £ (say 
for £ = 1 . . . 10) and the numerical evidence suggests that, for such small values of £, 
Q{£) is well-approximated by a formula of form Q(£) ~ A£ + B, i.e. 

Q(£)/£ ~A + B/£ (8) 

where A ~ 1.2 and B ~ 0.6. If it could be proven that there exist exact values for A 
and B such that this approximate formula holds for all I with an error term which tends 
to zero as £ tends to infinity, then this would of course be a stronger result than our 
proposition above and would, in particular, tell us that our value for ci, e' 1 / 3 ~ 0.71, 
can be improved to 1.2. However, we have been unable to prove this. 

2. Proof of Proposition 

We have already proven the upper bound above, so it remains to prove the lower bound. 
For this proof, we find it useful to write Q(£) as the product 

n 

Q(£) = Ylui(n) 

i=l 

% As far as we are aware, the notion of closest relevance to this question in the mathematical literature 
is the notion |2j of the transfinite diameter of a closed bounded region (which one shows, see again 
is equivalent to another notion known as the capacity of that region) of the (complex) plane. This is 
defined to be the supremum over all finite sets of points in the region of the l/n(n — 1) power of the 
product of all the distances between pairs of distinct points of the set. (One can show, for example, 
that the transfinite diameter of a disk is equal to its radius.) If we allow ourselves to generalize this 
concept of transfinite diameter to closed bounded regions of i/iree-dimensional Euclidean space, then 
our Q(£) may be seen to resemble one of the terms in the supremum which would define the transfinite 
diameter of a cube of side 2£+l (except that the 1/n 2 power is taken, instead of the 1 / n(n — 1) power) . 
However, our interest is in Q{£) itself and no supremum is to be taken. 
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where 

= Yl(\^-^\)^ 

where the product is over all j from 1 to n except that the factor with j — % is omitted. 
We shall also assume that the numbering of the vertices from % — 1 to n is such that 
the centre vertex (with coordinates, as introduced above, (£,£,£)) is numbered i = 1. 

Lemma. u)\ < ujj {j — 1, . . . , n). 

Proof. In terms of our coordinatization of our cube, we need to show that the product 
Rjuj 2 ,h) = 1 1 I Ui, h, h) - (M2, k) I 

(*l,*2,«3)^0'l ,32,33) 

over all distances from an arbitrary lattice point (ji, J2, J3) to all the other lattice points 
in our cube is greater than or equal to the product 

A^+M+M+i) = JJ \(£ + l,£ + l,£+l) - (i 1 ,i 2 ,i 3 )\ 

over all distances from the centre point, with coordinates (£ + !,£ + !,£ + I), to all the 
other lattice points in our cube. (To explain the notation here, Uj = (/^ij 2 ,.j3)) 1//n2 if 
(ji)j2, 33) are the coordinates of the point numbered j. The inequality of the lemma 
will then follow by taking the 1/n 2 power of each side of the above inequality.) We 
will prove this by exhibiting, for each j 2 , J3), a one-one onto mapping f(j u j 2 ,j3) fr° m 
(coordinatized) lattice points of our cube to themselves with the properties 

fUu32, j3 )((jlj2,js)) = (€ + M + M+ 1) (9) 

and 

l/0ij 2 j3)((M2,z 3 )) - (£ + l,£ + l,£ + l)\ < \(ii,i 2 ,i 3 ) - (ji,h,h)\ (10) 
In fact, given such a mapping we will immediately have 

Hh,j2,33) = 1 1 I U1J2, 3s) - (ii, h, is) I 

(h,i2,i3)^(jl,32,j3) 

> ) [ I f(h ,32,33) ( Ui , h ,33))- f(h ,32,33) ( (*i , *2 ,13) ) I 

{h,i2,i3)¥=Ul,32,j3) 

n |(€ + 1,£ + 1,£ + 1) - f U u32,33)((ii,i2, <s))| 

(u,i2,i3)^0'i J2J3) 

J] l^+l^+l^+l)-^,^,^)! 

= A^+i.m^+i) 

(11) 

where, in the penultimate equality, we have defined (^1,^2,^3) to be /(j 1 ,j 2 ,j3)((^i? ^2, h))- 
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It therefore remains to exhibit a mapping, /(ji j 2 3 - 3 ) , with the above properties. To 
do this, we define 

f(h,h,h)((h, 12, i 3 )) = + 1 - ii + M + 1 - J2 + h,t + 1 - J3 + z 3 ) (mod 2£ + 1) 

where we use a non- conventional notion of addition modulo q in which q (mod q) is 
deemed to be q rather than 0. One may understand this definition geometrically as a 
rigid translation of all the points of our cube where, however, if the would-be destination 
of a point is outside of our cube, the point instead gets mapped to the counterpart- 
position in our cubic lattice to the position it would arrive at in a neighbouring cube, 
were our cube to be surrounded by similar neighbours in a cubic lattice arrangement. 

Clearly this is a bijection and satisfies (jHJ). To show it satisfies (fTOj) we calculate, 
using the obvious properties of our notion of "mod" : 

l/«xJ^-s)((*i» *2, is)) ~(£ + l,£+l,£ + l)\ = ([((£ + + h)(mod 2£ + l))-l- l] 2 

+ [((£ + l-j 2 + z 2 )(mod2£+l))-£-l] 2 
+ [{(£ + l-j 3 + i 3 ) (mod 2£ + l))-£- I} 2 ) 1 ' 2 
< ((Ji - H) 2 + {32 - i2? + (3s ~ * 3 ) 2 ) 1/2 = |(jr, 32, 3s) ~ (ii,t 2 ,<s)|. 

□ 

Before proceeding with the proofs of the two bounds in our proposition, it will be 
useful to define, for a given cube (i.e. cubic lattice of points) of side 2£ + 1, certain 
special sets of lattice points. First, we define the "fcth centre shell" (k < £) to be the 
set of lattice points on the surface of the cube of side 2k + 1 centred on the centre of 
our given cube (so that a cube of side 2£ + 1 would have a total of £ centre shells). 

Proposition (Lower Bound). For a cube of side 2£ + 1, we have 

9!A> e -\ weN (12) 

Proof. For such a cube we have that 
n = (2£+ l) 3 

Let us denote the number of lattice points in the kth centre shell by rife. This is given 
by 

n k = 6(2k + l) 2 - 12(2Jfe + 1) + 8 
We first use the lemma to estimate 



Q(^) = na; i (7i)>a; 1 (n) n 



i=l 



We can rewrite u)\ (n) by a relabelling of the lattice points 

i 

nixi- Xl i = ( nn i x - ) ( i3 ) 
i=i / \fe=i d=i 
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where x^ fc is the coordinate-triple of the d th lattice point in the fcth centre shell. Each 
lattice point in the fcth centre shell is at least a distance k away from the centre lattice 
point so we have 

Q(£) > a(£)» (14) 

where 

e 

a(£) = Y[k nk (15) 

k=l 

From (|14j) . it can be seen that to prove (|12|) it will be sufficient to prove 

a{£) > (e—Bi) = ( e —B£j (16) 

This can be proved by induction on £. For all £ G N, let P{£) be the proposition that 
(Unj) is true. We have 

a(l) = 1 > e~s, 

so -P(l) is true. So, if we can show P(£) =>• P(£ + 1) (W e N) then the proposition will 
have been proved to be true. From (|T5|) and fTH|) 

a(£ + 1) = a{£){£ + l) ne+1 > (e"^V (£ + l)"^ 1 

Thus to now prove our proposition it will be sufficient to prove that 

(2M-3) 3 



i \ n / i \ ( 

e~B£\ (f+lf'« > r e -3(^ + i)J 



This is easily seen to be equivalent to 

In ,( 132^+72^-1 \ 

e-<(l + ? ) 

However the above equation clearly holds so we have proved our inductive proposition 
(and hence the lower bound in our main proposition). 

□ 
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